Temporally multiplexed superposition states of continuous variables 
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We study non- Gaussian states generated by two-photon subtraction from a cw squeezed light 
source. In a cw scheme one can subtract two photons from the source with a designated time 
separation and can genarate temporally multiplexed superposition states of continuous variables. 
We numerically study the properties of these states in the light of bosonic interference in the time 
domain. In an appropriate temporal mode amplified kittens are produced in a region where the 
time separation is comparable with the correlation time of squeezed packets. 

PACS numbers: 



I. INTRODUCTION 

In optical quantum information processing there have 
been two kinds of approaches with discrete- variable (DV) 
and continuous- variable (CV) schemes. The former usu- 
ally consists of entangled photon states and photon coun- 
ters, while the latter consists of squeezed states and ho- 
modyne detectors. Each of these schemes is, however, 
only a part of control with full potential aspects of quan- 
tum field of light. The scheme combining these two will 
surely be promissing for exploiting universal quantum 
information processing for computing, communications, 
and metrology. 

Recent progress in this direction has been made as 
the generation of quantum superposition of mesoscopi- 
cally distinguishable states using a squeezed state, pho- 
ton counting, and homodyne detection [l|, |2|, |3|, y, |5|. 
These experiments are based on the photon subtraction 
from a squeezed state, as proposed by Dakna et al. [6|. 
In the scheme a small fraction of squeezed light is tapped 
via a beam splitter, and is guided into photon counters. 
The remaining light beam conditioned on photon clicks 
turns to be a highly nonclassical and non-Gaussian state, 
which was transformed from a squeezed state (Gaussian 
state) via a strong nonlinear process of photon counting. 
In this way quantum superposition of mesoscopically dis- 
tinguishable states (Schrodinger kitten states) were gen- 
erated by the single-photon subtraction. 

Next step is to produce two or more kitten states as 
well as to increase the size of cat states for realizing 
quantum information processing and quantum metrology 
05 S @] • If larger number of photons could be subtracted 
from a squeezed light, larger size of quantum superpo- 
sitions (Schrodinger cat states) could be produced, al- 
though experiments get more challenging [10]. Another 
method is to produce a larger cat state from two smaller 
cat states as proposed in [11]. A new method to produce 
larger cat states was recently demonstrated in laboratory, 
by using_ two-photon Fock state with conditional homo- 
dyning 
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In this paper, we study another avenue in this direc- 
tion, which is the generation of two kitten states with a 
designated time separation in the temporal domain by us- 
ing the two-photon subtraction from a cw squeezed light. 
Two kitten states are generated in a single spatial mode. 
They are of temporally two- mode, and generally entan- 
gled, referred to as temporally multiplexed kitten states. 
We numerically study the properties of these states in 
the light of bosonic interference in the time domain, and 
show that in an appropriate temporal mode amplified 
kittens are produced in a region where the time separa- 
tion is comparable with the correlation time of squeezed 
packets. The paper is organized as follows. In Sec. Ullwe 
introduce basic notions for our cw scheme, and provide a 
rough sketch of physics with an idealized lossless model. 
In Sec. mil we present a precise modeling in a practi- 
cal setting, show numerical results, and discuss bosonic 
interference occuring in the cross over region where the 
time separation is comparable with the correlation time 
of squeezed packets. Section HVl concludes the paper. 



II. PHOTON SUBTRACTION IN A CW 
SCHEME: ROUGH SKETCH OF PHYSICS 

The photon subtraction scheme is depicted in Fig. [H 
A small fraction of the squeezed beam in path A is tapped 
at a beam splitter (BS) with reflectance R^ guided into 
two photodetectors in paths B and C through an opti- 
cal filter and another beam splitter with reflectance i?i, 
and is then used as trigger signals for conditional photon 
subtraction. 

We consider the scheme with a cw squeezed light gener- 
ated by an optical parametric oscillator (OP O). as stud- 
ied extensively by Nielsen and M0lmer [1^, [ij, [ill • A 
typical feature of this scheme is that the trigger photode- 
tection takes place in much shorter time scale (<1 ns for 
Si APD) compared with the OPO time scale (~100ns), 
and hence interesting multimode kitten states can be gen- 
erated in the time domain [2|, 0] • This is in sharp contrast 
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FIG. 1: (Color online) Scheme of the two-photon subtrac- 
tion from a squeezed state. BSs are beam splitters with the 
reflectances of R and Ri . 



to the short pulse scheme where the photon counting time 
scale is longer than the pulse width [l|, 0, [5|, [l2| . 

The annihilation field operator in a cw scheme is de- 
noted by 



1 f^ 

d{t) = — / dna{iJo + ^)e-^^^°+^^^ 



(1) 



where coq is the center angular frequency of the spectrum 
of a laser source. The operator obeys the continuum 
commutation relation 

[d{uJo + ^), d\uJo + ^0] = 27rJ(^ - ^0- (2) 

The time dependent field operator d{t) is defined in the 
interval (—00,00), and obeys the commutation relation 



d{t),d^{t')] =S{t-t'). 



(3) 



Squeezed light fields can be conveniently described in the 
rotating frame about the center frequency cjq, 



1 r^ 

A{t) = a(t)e^^°' = — / dnA{n)e 
27r J_oo 



-iQt 



(4) 



where A{Q) = d{u;o + Q). 

Now let us roughly sketch an essential physics by as- 
suming a nearly lossless OPO squeezing described by the 
Bogolubov transformation 

s\A{n)SA = ^{n)A{Q) + u{n)A\-n). (5) 

The cw beam described by a pure state Sa\^a) is the 
input squeezed state for photon subtraction. 



A. Odd-number Schrodinger kitten by 
single-photon subtraction 

Let us first consider the single-photon subtraction with 
a trigger signal detected at t = ti in the detector in path 
B {Ri = in Fig. [1]). We assume that this detection 



takes place instantaneously compared with the OPO time 
scale, and the trigger beam is projected onto the vacua 
for other times. Denoting the BS operation with the 
reflectance R as Vab, the conditional (not normalized) 
state is then described by 



.(1) 



p'^i^} (X {OB\B{ti)VABSA\OAB) 



-VRA{ti)SA\OA) 

POO 

VRSa / dtA\t)iy{t - h) \0a) (6) 



where we have used the relation ([5|), and introduced 

1 f^ 

u{t -ti) = — / dne'^^'-''^u{n). (7) 

27r J_oo 

We have also assumed that R is small so that the renor- 
malization of the effective squeezing on beam A due to 
optical loss is negligible. 

This quantity represents the temporal correlation of 
squeezed photons. It is roughly given by 



where 



Kt) 



m 



-Mt) 



Coe' 



-Col*l 



(8) 



(9) 



is a normalized temporal mode function p^, [IJ, Il5| . Here 
Co/tt corresponds to the OPO resonant bandwidth, and 
e is the nonlinear coefficient of the OPO, propotional to 
the x^^^ interaction coefficient and the pump field ampli- 
tude. Defining the field operator localized in this tempo- 
ral mode by 



/CO 
dtA\t)ip{t - h) 
-CX) 



(10) 



the single-photon subtracted state can be written as 



Pew 



oc 



VRe 



SaA\ \0a) . 



(11) 



This represents a cw beam consisting of the squeezed 
single-photon state [li,JJ,] localized in the mode il){t — ti) 
and the squeezed vacua in the other time domains. 

The beam is filtered through the mode function '^(t — 
ti), and then sampled by homodyne detection [13]. In the 
homodyne output one can see that two distinguishable 
waveforms with phase difference of tt coexist as demon- 
strated experimentally in [l|, Q, Q. The constructed 
Wigner functions have negative part, showing the non- 
classicality of the states. This is a quantum superposition 
of mesoscopically distinguishable states as analyzed by 
Dakna et al. [6], the so called Schrodinger kitten states. 
In this case it is an odd-number kittens. 
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(12) 



One can immediately see the two extremes; 

(i) A (= 1^2 — ^i|) <C Co""^ where the two-photon sub- 
traction takes place in the single mode, producing an 
even- number kitten as originally proposed in ^6|] . 

(ii) A :^ (q^ where two odd- number kittens are gen- 
erated in separated packets '^(t — ti) and '^(t — ^2). Each 
kitten is the one generated in [2|, [3] • 

In the imtermediate region there is a cross over be- 
tween (i) and (ii), and both even- and odd-kittens exsist, 
being entangled over the two modes. This region is of 
our interest. 

Since the field operators A[ and Al have generally a 
finite mode overlap for finite A 

/a = {0a\A2A\\0a) 

dttl;{t-ti)tl;{t-t2) 

) 

= (l + CoA)e-«°^ (13) 

one needs an orthonormal mode set to describe the state. 

1. Unbiased modes 

We first consider an orthonormal mode function pair 
{^[/i(t), ^f/2(^)} which suits to analyze the region CoA > 
1, and defined by 



^{t-t2) = VP^^Ul{t)^V^-PA'^U2{t) (14) 



where 
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B. Entangled kittens by tw^o-photon subtraction ' ' 



In the two-photon subtraction, we set Ri = 0.5 (i.e. 
50%) in Fig. [TJ Suppose that a first and second photon 
are detected at t = ti in path B, and at t = ^2 in path C, 
respectively. The conditional state can be represented as 
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FIG. 2: Temporal shape of the unbiased modes (dashed lines) 
in the case of Co = 30 MHz. The solid lines are '0(t — ti) and 

V^(t-t2). 
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FIG. 3: Vector representation of the unbiased modes. 



We call them the unbiased modes because ^f/i(t) and 
^c/2(^) equally overlap with il){t — ti) and '^(t — ^2), re- 
spectively. Their temporal shapes are shown in Fig. [2] 
as well as the mode overlap configuration in vector rep- 
resentation in Fig. [3l The corresponding field operators 
are defined by 



satisfying 



dtA(t)'^Ui{t) 



[Aui,A^UJ]=^^^ 



(16) 



(17) 



The operators A\ are represented as 



(18) 



The two-photon subtracted state (not norahzed) is then 
represented as 



A{t2)A{ti)SA \0a) 



Sa 



;^ |1,1)+/a 
Co 



|2,0) + |0,2) 
V2 



(15) 



-ee 



-CoA 



0,0) 



(19) 



The ket vectors in the above equation are understood as 
|- • • , • • •) = |- • •)jj-^ \- • •)u2' Thus the state in terms of 
the unbiased modes is invariant under the mode permuta- 
tion. The A-dependence of the above state is completely 
governed by the coefficients in the bracket. Note that 
Sa is independent of A. The index A means the other 
remaining modes than those spanned by the Ul and the 
U2. The modes A will be filtered out to observe the 
two-photon subtracted states. 

The two-mode state is essentially an entangled state 
of the two kinds of kittens (the odd-number and even- 
number kittens), and the squeezed vacuum. In fact the 
two-photon subtracted state can also be represented as 
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(2) \ _ A{t2)A{t^)SA\QA) 

Pew , 
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|0i).(20) 



where 



\<Pe 



f(l + CoA)|2) + i=|0) 
so v2 



AT 



so 



,e-2CoA 2 



-1/2 



|(l+C.AP + i 



(21) 



The first term includes two odd-number kittens in the 
similar sense of Eq. (pT]) . As for the second term, 
when |0e) is squeezed, one has an even- number kitten as 
demonstrated by Dakna et al. [6]. (Note that this is not 
precisely a coherent-state superposition but a superposi- 
tion of other kinds of two mesoscopically distinct compo- 
nents with even- number photons. The detailed features 
will be presented in the next section numerically.) 

The characteristic features are summarized as follows: 

(i) The first term in Eq. ([2Q|) includes the two single- 
photon-state kernels, each comes from the distinguishable 
photon subtraction from each squeezed light packet in 
i^{t — ti) and ilj{t — ^2). This corresponds to the two 
odd- number kittens. (See Eq. (pT]) .) 

(ii) The second and third terms with \(j)e) in Eq. ([2Q|) 
include the two even- number-state kernels. They come 
from the indistinguishable photon subtraction from the 
overlapping components between the two localized pack- 
ets in '0(t — ti) and '^(t — ^2), and their weights decreases 
as A. The |^e) being squeezed, the resulting state is an 
even-number kitten. 



FIG. 4: Temporal shape of the biased modes (one-dotted 
lines) in the case of Co = 30 MHz. The ^+(t) is symmetric, 
while the ^-(t) is asymmetric. The solid lines are ip(t — ti) 
and ip(t — t2). 
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FIG. 5: Vector representation of the unbiased modes. 



2. Biased modes 

Let us then consider the other important mode func- 
tion pair, which suits to analyze the region Co^ < 1, and 
defined by 



*±w = 



^(t-t2)±^(t-ti) 

\/2(1±/a) 



(22) 



Their temporal shapes are shown in Fig. [4] as well as the 
mode overlap configuration in vector representation in 
Fig. [5]. The ^+(t) is the mode that extracts the features 
due to the overlap of the packets '^(t — ti) and ip{t — ^2), 
while the ^-(t) is the mode that extracts the differential 
features of these localized packets. We refer them to 
as the biased modes. Defining the field operators in the 
biased modes by 



/CO 
dtA{t)^±{t) 
-00 



(23) 



|0+) the two- 



and the number states by |n+) = 

photon subtracted state can either be represented as 

A{t2)A{h)SA \0a) 



Sa 



e^ fl + lA 



CoV v^ 



2+,0_ 



1-/A 

V2 



0+,2_ 



+ee' 



-CoA 



o+,o_ 



lo^)- 



(24) 



Thus in the biased mode the output state is a squeezed 
entangled- state of 2-photon and the vacuum states over 
modes ^±(t). This is actuahy an entangled state of even- 
number kitten state and the squeezed vacuum state. For 
CoA<l 

A{t2)A{h)SA\0A) ^ ^A+fV2^|2+) + e|0+) 

V Co 

0Sa-\O-)^S^\O^) (25) 

where we have approximately decompose the squeezing 
operator Sa into a product of the operators for the modes 
^+(t), ^-(t) and others. Most of photons are contained 
in mode ^+(t) while mode ^-{t) is nearly the vacuum 
state. Actually the squeezing in the aymmetric mode 
^_(t) is very small. This feature holds for (qA ^ 1, and 
for certain finite A's, an amplified even-number cat state 
is produced in mode ^+(t) as shown in the next section. 
For (^0 A ^ 1 on the other hand photons are distributed 
equally to modes ^+(t) and ^-(t) as 



A{t2)A{h)SA \0a) 

-^ ^=— ^A+ ^ Sa- 

v^Co 



2+,0_)-|0+,2_) 



(26) 



Thus the two-photon subtraction from a cw squeezed 
beam generates temporally multiplexed Schrodinger kit- 
tens over the two temporal modes. In the next section 
we show the features described above by calculating the 
Wigner functions based on the precise modeling in a prac- 
tical setting. A new formalism to treat only the two 
modes of interest with the effective super-operators and 
detailed analysis of the states are presented elswhere 
119|. 



III. PHOTON SUBTRACTION IN A CW 
SCHEME: PRECISE MODELING 

A. Model 

In this section we extend the idealized lossless model 
described above to a practical case, and provide the for- 
mulae and numerical results for future experiments. In a 
practical setting of OPO, the pure-state argument in the 
previous section cannot apply because the two tempo- 
ral modes have generally quantum correlations with the 
other modes (other neighboring squeezed packets). This 
makes the relevant states mixed more or less. One needs 
to calculate density operators, whose details are given in 
Appendix. 



A cw squeezed state generated by an OPO, p)^\ is 
fully characterized in the time domain by the correlation 
functions [13] 



Tr [p^2^A{t)A{t' 




where 



C(e) = Co + e, Co 



7r + 7L 



(29) 



with the leakage rates jt of the output coupler and jl 
of the cavity loss. The Co determines the resonant band- 
width of the cavity. In a typical experiment with a low 
loss x^^^ -crystal H, jt ^ 57MHz, jl ^ 1.2MHz, and 
Co/tt ~ 9.3 MHz (the full width at half maximum). 

A small fraction is taken from the squeezed state and 
is split into paths B and C as the trigger beams, resulting 
the three-beam squeezed state 



PABC = VbcVabPa^vIj^V^c- 



(30) 



The trigger beams are then measured by two on/off pho- 
todetectors with the time resolution T (~ Ins). Its 
mathematical modeling is described in [18]. Denoting 
the "on" -signal elements of the positive operator- valued 
measure (POVM) of the on/off detectors as Ub and He, 
the output state in path A conditioned by the "on" sig- 
nals is given by 



1 



PA 



Pdet 



Tr^C [pABC^B ®^C 



where 



det 



TrABC {pABC^B (^IIc 



(31) 



(32) 



is the detection probability. 

This state is still a cw beam, consisting of multimodes. 
We are, however, interested in the two modes; the un- 
biased modes {^f/i(t), ^f/2(t)} and the biased modes 
{^+(t),^_(t)}. Denoting these sets as {^kif)} we in- 
troduce the annihilation operators 



dtA(t)^^k(t). 



(33) 



Let the quadrature amplitude and phase operators be 

It 5^^^!^. (34) 



f Ak + Al 



V2 



V2i 



In the homodyne channel one can observe a reduced state 
from pA into a desired mode by filtering the homodyne 
current into 



Xk{<P) 



Ake-'^ + Ale'^ 



V2 



(35) 




FIG. 6: Squeezing and antisqueezing characteristics of the 
tempral mode ip{t) = VT^ ■^'*' ^^ ^ function of e/Co- The 
sohd (dashed) hues are for / = 10, 30, and 50 MHz, from the 
bottom (top). 



and by constructing the Wigner functions with respect 
to this quadrature ampHtude. 



B. Numerical results 

In Fig. [6] the squeezing and antisqueezing characteris- 
tics of the input beam are shown for the tempral mode 
V;(t) = x/Je--^!^! for / = 10, 30, and 50 MHz. Since 
the squeezing itself is larger at lower frequencies, filtering 
with wider tempral shape observes larger squeezing. But 
the photon-subtracted squeezed state in the cw scheme 
is most appropriately observed for a characteristic fre- 
quency of the OPO / - Co -^ 30MHz (the thick solid 
and dashed lines in Fig. [6]). The reason for it includes a 
fact that the photon counting by Si APD takes place in 
a time scale much shorter than the squeezing correlation 
time. 

Figure [71 shows the Wigner functions (see Appendix for 
the formulae) and their contour plots; a is for the biased 
mode ^+(t), b is for the biased mode ^-(t), and c is for 
the unbiased mode ^f/i(t). The nonlinear coefficient is 
taken as e/Co = 0.3, corresponding to -3.6 dB squeezing 
for the localized mode '^(t) = vTo^~''°'^' • "^^^ reflectance 
of the tapping BS is set to R = 0.05 (i.e. 5%). The 
quantum efficiency of the homodyne detector is taken as 
r]H = 0.96 in Eq. (|A38p . For the on/off detectors, the 
overall detection efficiency of 77 = 0.6 and the fake trigger 
rate of u = 10~^ are assumed in Eq. (|A23p . Since the 
Wigner functions are identical for both of the unbiased 
modes, only the ones for ^f/i(t) are shown. From left 
to right A = 1, 30, 50, and 250 ns. Here note that the 
original two-mode state is essentially an entangled state. 
So a reduced state to a particular mode is more or less 
a mixed state. But one can see typical features of kitten 
components behind. 

The left two columns are typical cases of Co^ ^ 1- Toi" 
A = 1 ns an even-number kitten is generated in mode 



^+(t), while a nearly vacuum state is in mode ^-(t). 
The reduced state to mode ^f/i(t) seems close to a ther- 
malized squeezed vacuum. For A = 30 ns an ampli- 
fied even-number kitten appears in mode ^+(t) with a 
deeper negative Wigner function. The reduced state to 
mode ^f/i(t) becomes close to a mixture of two Gaussian 
states. 

The right two columns are typical cases of Co A > 1- 
For A = 50 ns the reduced state to mode ^+(t) still has 
a negative Wigner function, and larger amplitude. The 
state seen in mode ^f/i(t) shows a dip at the phase space 
origin, which is a typical feature of odd-number kitten. 
Actually VF(0, 0) directly reflects the parity of dominant 
photon number, namely a dip for odd and peak for even. 
For A = 250 ns one has two odd- number kittens as seen 
in mode ^[/i(t). In terms of the biased modes ^±(t) 
the states reduces to identical mixed states of the even- 
number cat and the squeezed vacuum. 

Figure [8] shows the case of A = 65 ns. The upper pan- 
els a and b show the Wigner functions and their contour 
plots of the states reduced to the unbiased modes, while 
the lower panels c and d show the ones of the states re- 
duced to the biased modes. For this time separation the 
reduced states in both kinds of modes show the nega- 
tive Wigner functions. The panels a and b correspond to 
the odd- number kittens, while the pannel c and d corre- 
spond to the squeezed vacuum and the even- number cat. 
(These reduced states are more or less mixed ones due to 
the two- mode entanglement behind.) Thus by changing 
the filtering functions one can access different features of 
the components. One possible way of viewing the change 
from a and b to c and d in Fig. [8] is a constructive bosonic 
interference of the mesoscopic superposition states in the 
sense that odd-number-rich components seen in modes 
^f/i(t) and ^f/2(^) condensate into mode ^+(t) construc- 
tively, producing an amplified even- number-rich cat. On 
the other hand, in mode ^-(t) the destructive interfer- 
ence occurs, leaving the vacuum-like state. 

In Fig. Owe depict average photon numbers as a func- 
tion of A for five kinds of squeezing levels, from the bot- 
tom e/Co = 0.1 (-1.2dB), 0.2 (-2.4dB), 0.3 (-3.6dB), 0.4 (- 
4.7dB) and 0.5 (-5.7dB). See Fig.[6]for the corresponding 
squeezing levels for / = 30 MHz. (The formulae are given 
in the last part of Appendix.) The left is for the biased 
modes ^+(t) (solid line) and ^-(t) (dashed line). The 
right is for one of the unbiased modes. Since (ni) = (712) 
for the unbiased modes, the right panel directly show how 
the total average photon number increases and saturates 
to a certain level as (qA increases. 

An interesting feature is the peak structure in mode 
^+(t) seen in the left panel. Most of photons accumu- 
late in mode ^+(t) around CoA ~ 1 ~ 3, depending on 
the squeezing levels. On the other hand, the average 
photon number in mode ^-(t) is suppressed in this time 
range. Roughly speaking amplified even- number kittens 
appear in this time range. A typical Wigner function 
for larger squeezing e/Co = 0.5 with A = 43 ns is shown 
in Fig. [TOl This is close to a quantum superposition of 
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FIG. 7: Wigner functions and their contour plots, a for the biased mode ^+(t). b for the biased mode ^-(t). c for the 
unbiased mode ^c/i(t). From left to right A = 1, 30, 50, and 250ns. 
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FIG. 8: Wigner functions and their contour plots for the two sets of modes in the case of A = 65 ns. 
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0.4 and 0.5. The left is for the biased modes ^+(t) (solid line) and ^-(t) (dashed line). The right is for one of the unbiased 
modes. 
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FIG. 10: Wigner function and its contour plot of the state 
in mode ^+(t) around the peak of (n+), corresponding to 
e/(o = 0.5 and A = 43 ns. 



distinct squeezed states with phase difference of tt. In- 
terestingly this is quite similar to the cat state generated 
by four-photon subtraction, which was predicted in [6|. 
Precisely speaking the maximum quantumness of the cats 
including the negativity of the Wigner function and the 
state purity does not necessarily coincide to the peak of 
(n+). The detailed analysis will be presented in another 
paper p^ . 

The process of producing the amplified kittens is noth- 
ing but the constructive bosonic interefence caused by 
nonlinear process of photon subtraction. This is directly 
seen as the superposition in Eq. ([M|) . The mechanism 
behind this consists of the following two aspects. The 
first one is the indistingushable two-photon subtraction 
from the temporally overlapped squeezed packets. This 
aspect is represented by the term 6'^A2A| \0a) in Eq. 
([12]), or equivalently the first and second terms in Eq. 
([M|) . Note that the two-photon subtraction leaves the 
state that can be represented as the squeezed two-photon 
states. The second one is the existence of the squeezed 
vacuum as the second term in Eq. ([12]). This corresponds 
tothetermS'|0+,0_) [=6'|0,0)] in Eq. ([24]) [Eq. (p^]. 
This term allows one to tune the coefficients of super- 
position between the squeezed two-photon state and the 
squeezed vacuum. The second aspect makes a sharp con- 
trast with a similar scheme with nondegenerate two-mode 
squeezed vacuum to generate two-photon state proposed 
in [l5[, where such a term does not exist. 

The state structure with the parameters A and the 
squeezing coefficient e explained above allows one to con- 
trol the bosonic interference in the time domain, and gen- 
erate variety of quantum superposition states. 



IV. SUMMARY 

We studied non- Gaussian states generated by two- 
photon subtraction from a cw squeezed light source. A 
trigger photon click specifies a certain temporally local- 
ized mode in the remaining squeezed beam. The two- 
photon subtraction from overlappping squeezed pack- 



ets generates temporally multiplexed Schrodinger kittens 
with a designated time separation A. They are gener- 
ally an entangled states of even- and odd-number kittens, 
and the squeezed vacuum over two temporal modes. We 
have numerically studied these states by calculating the 
Wigner fuctions of the reduced states into each mode of 
the two kinds of sets. In the biased mode ^+(t) ampli- 
fied kittens are produced in a region A ~ (^q~^ where the 
time separation is comparable with the correlation time 
of squeezed packets. This is due to constructive bosonic 
interference in the time domain. 

As for furure perspective, it is interesting to consider 
ways of combining temporally multiplexed cat states with 
spatial degree of freedom to implement resource-effective 
quantum computation, and to generate useful nonclassi- 
cal resources for metrology. It is also interesting to apply 
the similar operation with overlapping squeezed packets 
to other kinds of bosonic systems in other degrees of free- 
dom such as the spatial and/or momentum domain. 
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APPENDIX A: MATHEMATICAL FORMULAE 
FOR TWO-PHTON SUBTRACTION 

In this Appendix we present methematical formulae to 
analyze the two-photon subtracted squeezed states in the 
cw scheme. 

The trigger channel consists of on/off detectors with 
the the time resolution T (~ Ins). The trigger photon 
mode is assumed to be rectangular with the duration T 



1 



c^k{t) 







exp(- 



27rkt, 



T T 

2 - - 2' 

otherwise. 



(Al) 



Since ^o^ ^ 1 only the lowest mode (poit) is excited in 
the photodetectors, and the weights of the other higher 
modes are negligible [18]. The trigger fields at ti- and 
t2-segments are defined by 



B 



C 



VT Jti-T/2 
I nt2+T/2 

\IT Jt2-T/2 



dtB(t) 
dtC(t). 



(A2) 
(A3) 



Let IIb and lie be the "on" signal elements of the pos- 
itive operator- valued measure (POVM) of the on/off de- 
tectors given in [18.] . 
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The output state in path A conditioned by the "on" 
signals is given by 



1 



PA 



Pde 



TlBC {pABC^B ^n, 



C 



where 



det 



Tr ABC {pABC^B ® nc) 



(A4) 



(A5) 



is the detection probabihty. Here note that we discard 
the information of the "off' signals. 

This state is still a cw beam, consisting of multimodes. 
We are, however, interested in the two modes; the un- 
biased modes {^ui(t)-,^U2(t)} and the biased modes 
{^+(t),^_(t)}. Denoting these sets as {^k(t)} we in- 
troduce the annihilation operators 



dtA{t)^^k(t). 



(A6) 



Let the quadrature amplitude and phase operators be 



their vector notation be 



X2 



^/2^ 



Pi 



(A7) 



(A8) 



and for beam B and C be 


Xb = 


V2 ' 


Xc = 


c + ct 

/7T ' 



^ ' 



Pc = ^=-. (AlO) 



The three-beam squeezed state in the relevant modes 
is conveniently described by the characteristic function 

C{pABC] UA, VA, Ub, Vb.Uc, Vc) 



Tr ABC S PABcexp 



^(^uaXa+VaPa 



^ubXb + vbPb + ucXc + vcPc) 



exp 



--^ur(-e)u 
4 ^ . 



exp 



--Vr(e)v 
4 ^ . 



where 



ua 




_ / ui 



va 




- / ^1 



(All) 

(A12) 
(A13) 



The covariance matrix is given by 



AA 



r(e) 



(e) 



^r^^R'liB{^) 



^iic(e) 



VR'iici^) ^/{l-Ri)Rib+-{e) l + Rib++{e) 



(A14) 



where 

TAA{e)=Rl + {l-R) 



Mn(e) Mi2(e) 
M2i(e) M22(e) 



(A15) 



M.(e) = 4,-^£ 






dt 



di'*fe(t)*;(*')e"'^^'^'*"*'' 

) 

(A16) 






(A17) 



/CX) 



ef (e) = Vii-R)R^^^ 



(A18) 

) 

(A19) 



6..(e) = -i?^ 



-(e) 



!)++(e)e' 



-C(e)A 



(A20) 

(A21) 

The characteristic function of the two-photon sub- 
tracted squeezed state is given by 



C(/9a;ua,va) 
■ Tr ABC 



Pdet 



1 



dUB / dVB / <^'^C / <^'^c 



(2^?Pdet 

xC{pABC; UA, VA, ^iB, ^B, ^C, ^c) 

xC(nB; -^/B, -VB)C{nc; -uc, -vc). (A22) 

The characteristic function of the "on" signal POVM el- 
ement is given by 



C{IVb',ub,vb) = 2tt8{ub)5{vb) 



-exp 



At] 



■(ii| + 4)jA23) 
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where rj is the overall detection efficiency, and u is the 
fake trigger rate. Substituting it into Eq. ()A22p we have 



c;(/)a;ua,va) = 2^-f; — ^^p 



=1 



det 



X exp 



V^r(^)(e)v^i 



where 



AA(i) = 1 
A/-(2) = - 



(A24) 



(A25) 



2e- 



r/V^(2)(-e)5(2)(e) 

26-"^ 
r?v/&(^)(-e)6(3T(^ 
4e-2- 



^2^5(4)(_g);,(4)(g) 



(A26) 



1 

2 



6(1) (e) 

6(2)(e) = -i?i6++(e) 

V 

6(4)(g) = 6(3)(e)6(2)(g) _ (1 _ R^)R^bl_{e) (A30) 



(A27) 
(A28) 

(A29) 



det 



= T.^'' 



i=l 



.(0 



,-(i) 



r^V(e)=i?+(l-i?)M«(e) 



6W(e) 



(A31) 



(A32) 






(A33) 

Ri^^ieKfie) (A34) 

(l-i?i)4^(e)Cf(e) (A35) 

(l-i?i)5(2)(e)Cf(e)^f(e) 

2(l-i?i)i?i5+_(e)Cf(e)?F(e) (A36) 



In real experiment the output state suffers from prop- 
agation loss, and imperfect efficiency of homodyne detec- 
tion. Taking these factors into account by the effective 
transmission coefficient tjh the observed state p^ is char- 
acterized by 



C'(/5^;ua,va) 



2^7^ exp 



det 



X exp 



.(i) 



^u^r^^(-6)u^ 



'vAT%\e)vA 



(A37) 



where 



T%\e) = {l-rjH)I^mr^'He) 



(A38) 



The output state also suffers from invasion of external 
noise photons from various parts, which degrades quan- 
tum correlations. This effect can be modeled by a mix- 
ing of thermal vacuum. Mathematically this can be 
taken into account by replacing the factor e^r/Ci^) i^ 
Eqs. (|AT6l) . (|AT8l) . (|AT9l) . and (|X2Ql) with (e-e^)7T/C(e) 
where Cx represents the degree of external noise photons. 
The two-mode Wigner function is then given by 



= v^ 

Z^ -,r2j 



exp 



*xr^'(-e)-ix 



n^Pde 



r^^(-e) 



exp 



*pr^\6)-ip 



.« 



(A39) 



r^^(e)l 



It is hard to depict the two-mode Wigner function nu- 
merically in a graph. One way to understand the state 
properties is to see the Wigner function reduced to an 
appropriate single mode ^/e(t) {k = [/I, [/2, +, or — ). It 
is given by 



Wix,p) = 



4 

E 



^i^) 



exp 



p^ 



■ Hkk 



i-e) 



■ Hkk 



(e) 



.(0 



(A40) 
The photon number in each mode can be calculated by 
the characteristic function by 



- /it 



(hk) = (AlAk) 






dv 



^3 ]C{pA;uA,yrA) 



J ua=va=0 



4Pde 






rS(-e) + rSW 



(A41) 



This can be reduced to explicit expressions for small 
CoT < 1 as 



(n±) 



2(l±7A)[z2(l+/i) + e-2A] 



(A42) 



where z = Co/e 
g[^\A) 



(l + /i)5i^)+2(l±/A)/i*^' (A43) 
G^^\A) ^ V(A)^ffi^^+2(l±/Af 

±4(l±/A)V(A)/i='^' (A44) 



with 



A^' - ^(l±/.)±^e- 



.i"^ 



5(1±/A)±e-^(2A^ + — ). 



For the unbiased modes 

{nui) = {nu2) 



(A45) 
(A46) 
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-[(n+) + (n_)]. (A47) 
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